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Prospects for gravitational-wave observations of neutron-star tidal disruption in 

neutron-star/black-hole binaries 

Michele VallisncriB 

Theoretical Astrophysics 130-33, California Institute of Technology, Pasadena CA 91125. 
INFN, Sez. di Milano, Gruppo Collegato di Parma/ Universitd di Parma, 43100 Parma, Italy. 

(December 7, 1999) 

For an inspiraling neutron-star/black-hole binary (NS/BH), we estimate the gravity- wave frequency 
/td at the onset of NS tidal disruption. We model the NS as a tidally distorted, homogeneous, 
Newtonian ellipsoid on a circular, equatorial geodesic around a Kerr BH. We find that / t< j depends 
strongly on the NS radius R, and estimate that LIGO-II (ca. 2006-2008) might measure R to 15% 
precision at 140 Mpc (~ 1 event/yr under current estimates). This suggests that LIGO-II might 
extract valuable information about the NS equation of state from tidal-disruption waves. 

PACS numbers: 04.80.Nn, 04.30.Db, 97.60. Jd 



The equation of state of the bulk nuclear matter inside a 
neutron star (NS) is poorly understood JlJ. For example, 
candidate equations of state that are compatible with 
nuclear physics experiments and theory predict, for a 
1.4 Mq NS, a radius anywhere from about 8 km to 16 km 
f|. Thorne has conjectured that insights into the equa- 
tion of state might come from measurements of the grav- 
itational waveforms emitted by merging NS/NS binaries 
and/or tidally disrupting NS's in neutron-star/black-hole 
(NS/BH) binaries |,|. More recently, Newtonian mod- 
els of NS /NS mergers have given strong evidence that the 
merger waves do carry equation-of-state information, but 
for NS/NS are emitted at frequencies (~ 1400-2800 Hz) 
too high for measurement by LIGO-type gravity-wave 
interferometers Iffi. In this paper, we show that the 
prospects for NS/BH measurements are much brighter. 

Central to these prospects is the question of whether 
NS tidal-disruption waves lie in the band of good inter- 
ferometer sensitivity (for LIGO-II, ~ 30-1000 Hz 0; see 
Fig. [I]). Numerical modeling of NS tidal disruption in 
NS/BH binaries is only now getting underway [0 and 
has not yet included computations of the emitted gravity 
waves or even their frequency bands. As a result, the 
best frameworks now available for estimating the tidal- 
disruption gravity- wave band are highly simplified, quasi- 
analytic models by Shibata || and by Wiggins and Lai 
[10 1, which represent the inspiraling NS as an irrotational 
[11 1, incompressible or polytropic Newtonian ellipsoid, 
moving on a circular, equatorial geodesic orbit around 
a Kerr BH, and being tidally distorted by the Kerr Rie- 
mann tensor. For simplicity we focus on Shibata's ho- 
mogeneous models, and then appeal to the polytropic 
models for evidence that compressibility has only small 
effects. 

In Shibata's analysis, the NS gravitational field, its 
centrifugal potential, and the Newtonian tidal poten- 
tial constructed from the Kerr Ricmann tensor are all 
quadratic functions of position. As a result, a class of 
equilibrium solutions are the classic irrotational, homo- 
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FIG. 1. Plot of (square root of) noise spectral density for 
different LIGO configurations: (1) LIGO-I; (2) LIGO-II wide- 
band; (3), (4) LIGO-II narrowband centered on 500 and 850 
Hz. Curves (l)-(3) are from [Q; curve (4) was produced by 
K. A. Strain using the same detector specifications as in Q. 



geneous Roche-Riemann ellipsoids ]r2] ]. Given a choice 
of the binary parameters M, a and r (the BH mass and 
angular momentum per unit mass, and the orbital sepa- 
ration, i. e., the Boyer-Lindquist radius of the geodesic), 
there is a one-parameter family of such NS models with 
density p ranging downward through the family to a min- 
imum p cv (M, a, r). 

We model the inspiraling NS as one of Shibata's irrota- 
tional ellipsoids, identified by its mass m, and its density 
p or mean radius R = (3m/47rp) 1 ^ 3 . In our simple frame- 
work, the uncertainty in m(R) embodies the uncertainty 
about the NS equation of state. We describe the inspiral 
as a sequence of circular, equatorial Kerr geodesies that 
shrink inward until the NS reaches the innermost sta- 
ble circular orbit, r = 7"j sco , or begins to tidally disrupt 
[which happens at the radius r t d where the star's density 
p matches the critical density p cr (M, a,r t d)]- 

The Kerr geometry provides a one-to-one correspon- 
dence between the orbital radius rtd an d the gravity- wave 



1 



frequency ftd at which tidal disruption begins: 



ftd(M,a,r td ) 



1 



(1) 



(here and below we set G = c = 1). It is this ftd that 
LIGO-II can measure. Having measured ftd and deter- 
mined the masses M and m from the observed inspiral 
waveforms fDj] , one can compute rtd and then the NS 
density p — p cr (M,a,rtd) and the mean NS radius R. 
Thereby the LIGO-II observations can determine a point 
on the NS mass-radius curve m(R), which represents the 
NS equation of state in our simplified analysis. Even 
one such point could give valuable information about the 
real NS equation of state, and several such points could 
determine it remarkably well [ fTi] . 

To estimate the accuracy with which LIGO-II might 
determine the NS radius R, we need the explicit rela- 
tionship between R and the disruption-onset frequency 
ftd- More precisely, we need R(m, M, a, /td), which can 
be derived as follows: (i) r t d(M, a, / t d) is obtained by in- 
verting Eq. P); (ii) p CT (M, a,r t d) is obtained by solving 
Eq. (3.9) of |J] for the ratios of semiaxes of the equilib- 
rium configurations, and then extremizing Eq. (3.10) of 
||, in which O 2 = M/(irpr 3 ); (iii) then R is obtained as 
R = [3m/47T/9 cr (M, a, rtd)] 1 ^ 3 - The result has the form 



R(m, M, a, / td ) = m 1/3 M 2/3 D 



M 



/tdM 



(2) 



where D is a dimensionless function with remarkably 
weak dependence on a/M jl5|]. This R(ftd) is shown 
in Fig. H for various M, for a/M = 0.998 (the curves 
for other a/M are almost identical to these), and for 
m = 1.4 M Q @]. The radii shown, R = 8-16 km for 
m = 1.4 M©, correspond to the range of predictions by 
plausible NS equations of state (2) . The curves in Fig. g 
are well approximated by the formula (with G = c = 1) 



R 



m 1 / 3 M 2 / 3 



0.145 (/ td M)-°- 71 
0.069 (/tdM)- 95 



for ftdM < 0.045, 
for f td M > 0.045. 



(3) 



Although the BH spin parameter a has negligible influ- 
ence on the function i?(/ t d), it strongly influences the ra- 
dius rj sco of the innermost stable circular orbit [[l7j . If the 
NS is still intact when it reaches ri SCO , it then will plunge 
rapidly into the BH and the tidal-disruption waves, if any, 
are likely to be so weak and short-lived as to be useless for 
measuring NS properties. Thus, there is not much hope 
of measuring tidal disruption unless ftd < /plunge = [Eq. 
(|Tj) with r td replaced by r isco (M,a)]; i. e., unless / td is 
left of the relevant big dot in Fig. ||. 

Figure || and the above discussion show that (i) for a 
wide range of realistic parameters, tidal disruption oc- 
curs before the plunge begins, and (ii) for all realistic pa- 
rameters except a very narrow range (M < 10 Mq and 
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FIG. 2. NS radius R vs. disruption-onset frequency / t d, 
for m = 1.4 M Q and M = 2.5-80 M Q . The black dots, 
parametrized by a/M, specify the onset of plunge into the 
BH; tidal disruption is measurable only for / t( j left of the 
plunge point, i. e., for R above it. Negative a/M indicates 
retrograde NS orbits. 



R < 10 km), the tidal-disruption waves fall in the range 
of good LIGO sensitivity, / < 1000 Hz. The Lai- Wiggins 
polytropic NS models Jl0[ ] give similar curves and con- 
clusions: for polytropic indices n = 0.5 and 1.0, which 
approximate NS equations of state, the R(ftd) curves are 
displaced upward in frequency from those of Fig. ^| by a 
mere ~ 50 and ~ 100 Hz. 

Turn now to an estimate of the accuracy to which 
LIGO-II could measure / t d (and then R) using Wiener 
optimal filtering \ 18 jlT|| . The measured gravity- wave data 
stream g(t) is compared to a set of theoretical inspiral 
templates h{9 % \t), indexed by the parameters Q % of the 
binary; a "best fit" 9 l is found which maximizes the like- 
lihood of observing g(t) given a "true" signal h(0 l ; t), and 
given a statistical model of the detector noise [a Gaus- 
sian [^o| random process with zero mean and spectral 
density S n (f)]. For strong enough signals, 9 l will have 
a gaussian distribution centered around its "true value" 
9 l , with covariance matrix iflSfl 



di = (r- 1 ) 



p.. 



dh 



dh 



(4) 



where the "inner product" (...) is defined for any two real 
data streams g(t), h(t) in terms of their Fourier trans- 
forms <?(/), h(f) by 



(9,h) 



df 



~g(f)h*(f) 

Sn(\f\) ' 



(5) 



Because so little is known about the tidal disruption 
and our NS models are so crude, we use the simplest 
of templates in our analysis: slow-motion, quadrupo- 
lar waveforms for point particles in circular, Keplerian 
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orbits with quadrupole-governed inspiral. The Fourier- 
transformed waveform, squared and averaged over binary 
directions and orientations, is given by plj 

{lhbl2) = 30^ J^tfW 0(/plungc ~ J) ' (6) 

where \i and Mt are the reduced and total masses, d is 
the distance to the binary, and the step function shuts 
off the signal at the onset of plunge. 

For typical observations, optimal filtering of the inspi- 
ral signal should give good estimates of M and m [ [l3|| . 
We therefore assume that the accuracy in measuring R 
is limited only by the uncertainty of /td [p2f - The es- 
timation of /td depends heavily on the details of the 
tidal-disruption waveforms, which are largely unknown. 
However, it is reasonable to expect tidal disruption to be 
a sudden event that significantly weakens gravity-wave 
emission within a few dynamical time-scales after / t d has 
been reached |^3| . Correspondingly, we employ a toy 
model where the inspiral waveform of Eq. (|^) dies out 
over a frequency band (/ td , /td + Sf): 

( h b (f) if / < /td, 

hd(f) = I hifm^) if /td < / < /td + Sf, (7) 
I if / > /td + Sf, 

where 0(a;) = 1 — x (linear decay), or Q(x) = I0~ x 
(exponential decay). The standard deviation of the "best 
fit" / td is given by Eq. (§ as A/ td = [T /td/td (M]" 1 / 2 . 

We have evaluated A/ td numerically, using the signal 
model from Eqs. (||), (Q) and the inner product ^ with 
the LIGO-II noise curves S n (f) of Fig. |l|. We have then 
computed the la range of the NS radii R from the re- 
lation R± = R(m,M,aJ td T 2A/ td ) [Eq. (§)]. The un- 
certainty in R, defined as Ai? = (R + — i?._)/2, scales 
roughly linearly with d [because A/ td is proportional to 
d through Eqs. (||), @], and is quite sensitive to the 
choice of the shutoff model [it scales roughly as (Sf) 1 / 2 
and is lower for the exponential decay than for the lin- 
ear one] . In Table Q we report the fractional uncertainty 
AR/R, averaged over the range 10 km < R < 15 km, for 
choices of parameters motivated by the following. 

The NS mass m was set to be 1 .4 M Q |L6| . The distance 
d and the BH masses M were chosen to represent two dif- 
ferent scenarios: (i) low-mass BH's, with M = 2.5 M Q at 
65Mpc (~ one merger/yr according to Bethe and Brown 
Hi); (ii) higher mass BH's, with M = 10, 20, and 40 M Q , 
at 140 Mpc (massive main-sequence binaries are thought 
to produce NS/BH binaries with M ~ 10 M and coa- 
lescence rates up to ~ one event/yr out to 140 Mpc, but 
possibly much less |2j|; capture NS/BH binaries formed 
in globular clusters might have M as large as hundreds 
of Mq [^6), but with exceedingly uncertain rates). Fi- 
nally, we considered three different gravity-wave shutoff 
models: (i) an optimal-precision model with linear decay 



TABLE I. Fractional uncertainty AR/R, averaged over 
the range 10 km < R < 15 km. Rows: BH masses; columns: 
gravity-wave decay models and detector noise curves [labeled 
by (2)-(4) as in Fig. [|. No quote is given if AR/R > 25%. 

AR/R (%) for 10 km < R < 15 km 

lin. decay exp. decay e. d., n = 1 lin. decay 

*/ = /td/6 5/ = /td/2 5f = /td/2 5f = /td/2 

M\S n (2) (3) (4) (2) (3) (4) (2) (3) (4) (2) (3) (4) 

2.5 M Q a 12 18 c 8 20 - 13 24 c - 25 21 c - 17 

lOM b 14 17 c 10 23 17 25 c - 21 25 c 19 

2OM b 10 16 10 16 14 c 16 23 14 22 25 c 16 

40M Q b 7 6 11 11 10 19 20 10 c 20 17 19 23 

a At 65 Mpc. b At 140 Mpc. c For 12 km < R < 15 km. 

and Sf = /td/6 (the lower limit set by the uncertainty 
principle on the frequency spread of waves emitted dur- 
ing 3 orbital periods, supposedly a typical time-scale for 
complete disruption |p3| ); (ii) a fiducial model with ex- 
ponential decay and Sf = /td/2 (a scaling supported by 
numerical calculations of tidal-disruption waveforms in 
NS/NS binaries ||); this model was also used to evalu- 
ate errors for n — 1 polytropes; (iii) a conservative model 
with linear decay and Sf = /td/2. 

The estimates for our fiducial decay model suggest that 
R may be determined with a precision of ~ 15% using 
the 850 Hz-narrowband LIGO-II configuration [curve (4) 
of Fig. , and with a somewhat worse precision for wide- 
band LIGO-II [curve (2)]. If the optimal-precision decay 
model is correct, the error might be as low as ~ 6-10%. 
The usefulness of the 500 Hz-narrowband interferometer 
[curve (3) of Fig. [j] is limited to the the heavier BH's or 
to the larger NS's, which have lower / t d- Our estimates 
are inferior for the Lai- Wiggins compressible polytropes 
]l0| examined in the least favorable case (n = 1), and for 
the most conservative decay model; even then, an 850 Hz- 
narrowband LIGO-II might be able to provide significant 
information about R. 

The accuracy of our analysis is limited by several fac- 
tors. Sources of error in the frequency /td(^, M, a, R) at 
which tidal disruption begins to significantly change the 
inspiral waveforms include: (i) the use of the test-mass 
approximation for the NS orbit, when actually m M, 
especially for the low-mass Bethe-Brown case; (ii) the 
use of the Riemann tensor to compute tidal forces when 
the NS diameter is not, typically, small compared to the 
distance from the NS center to the horizon |g7| ; (iii) the 
idealization of the NS as a homogeneous or polytropic el- 
lipsoid; (iv) the fact that the point at which the observed 
waveforms show a clear deviation from a standard inspi- 
ral may actually come a few orbits earlier (due to tidal 
coupling) or later than / t d- 

Our method presupposes a reliable technique to distin- 
guish a plunge shutoff of the inspiral waves from a tidal- 
disruption shutoff. In fact, it seems likely that the tidal- 
disruption waveform will actually contain features that 
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not only distinguish it from a plunge shutoff, but that 
also carry equation-of-state information which is richer 
than in our crude model. For example, simulations || of 
tidal disruption in NS /NS binaries show a spectrum with 
an inspiral cutoff followed by a valley, a moderately sharp 
peak, and a cliff; however, the NS/BH case is likely to be 
different, and the issue will ultimately be settled only by 
detailed numerical simulations. 

Given these large uncertainties, our results can only 
be rough indications of the prospects for learning about 
NS's from tidal-disruption waveforms. They do, however, 
suggest that observations of tidal disruption in NS/BH 
binaries might be possible in ~ 2006-2008 with L1GO-11, 
and might yield useful insights into the NS equation of 
state. The success of this endeavor will require the devel- 
opment of better theoretical and numerical techniques for 
modeling NS tidal disruption and computing the depen- 
dence of the disruption waveforms on the NS equation of 
state; we strongly advocate such an effort. 
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K. Alvi. This research was supported in part by NSF 
Grants AST-973f698 and PHY-9900776. 
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